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Abstract 

The characteristic multi-dimensional integrals that represent physical quantities in random-matrix mod- 
els, when calculated within the supersymmetry method, can be related to a class of integrals introduced 
in the context of two-dimensional conformal field theories by Dotsenko & Fateev. Known results on these 
Dotsenko-Fateev integrals provide a means by which to perform explicit calculations (otherwise difficult) 
in random-matrix theory. We illustrate this by (i) an evaluation of the mean squared S-matrix elements 
for the Gaussian orthogonal ensemble coupled with M external channels, and (ii) a direct derivation of 
the asymptotic behaviour of the dynamical density-density correlator in the limit of large spatial and 
temporal separation for the Calogero-Sutherland model which, at certain couplings, is known to map 
onto the parameter-dependent random matrix ensembles 
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1 The GOE Problem with External Channels 



1.1 B ackgr ound 

Various problems in quantum physics involving classical chaos or disorder can be modelled by Hamilto- 
nians belonging to one of the standard ensembles of Gaussian-distributed random matrices Q coupled to 
(non-chaotic) external channels. The random Hamiltonian models the chaotic dynamics of the internal 
system — which may be the compound nucleus in a nuclear scattering problem a ballistic electron- 
microstructure in the shape of a classically chaotic cavity ^ , or a disordered mesoscopic wire exhibiting 
diffusive electron motion Q. In the example from nuclear physics, the external channels represent the 
open decay modes of the compound nucleus; while in the mesoscopic conductance problems, they enumer- 
ate the transverse modes of the electron wave-function which enter or leave the chaotic system through 
attached ideally-conducting leads. 

The central physical observable of relevance in such systems is the S-matrix, for which the external 
channels constitute the asymptotic states. In nuclear-physics applications, the cross-section is related to 
the square of S-matrix elements. In applications from condensed-matter physics, it is the conductance 
that can be obtained by summing the squares of the S-matrix elements describing transmission. This 
relationship arises from the Landauer formula, viz., 

M/2 M 

E {\Sab\' + \S,a\'}, (1.1) 
a=l b=M/2+l 

where the index a sums over the M/2 incoming channels and b sums over the M/2 outgoing channels. 
Here, g here denotes the dimensionless conductance, and is related to the physical conductance G by 
g = h/e^-G. There is, in fact, a close analogy between the statistical behaviour of cross-section and con- 
ductance for such random-matrix problems. For example, there appear parallel universal characteristics 
in the two types of system: Ericsson fluctuations in nuclear cross-sections |^ and universal conductance 
fluctuations in mesoscopic physics |^, ^ . These comparisons are discussed more fully in Ref . ^ 

Thus, in order to study the statistical properties of the cross-section or conductance, one is interested 
in calculating the ensemble averages of moments of S-matrix elements over the distributions of the random 
Hamiltonian matrices. 

The coupling of the random system with the external channels c = 1, 2, . . . , is fully characterized 
by the so-called sticking probabilities Tc, which range between zero and unity. The value of Tc can be 
loosely viewed as an impedance mismatch for a given channel, which interpolates between total instanta- 
neous reflection and perfect transmission for values of zero and unity, respectively. The Tc are derived as 
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expressions involving combinations of all the microscopic parameters of the lead-system interface. These 
quantities are generally unknown and one does not try to specify them. Rather, assumptions are made 
directly about the sticking probabilities, based on the macroscopic characteristics of the experimental 
set-up. In most cases, the channels are physically equivalent, and so it is usual to take all the sticking 
probabilities to be equal. Also, in mesoscopic conductance problems, in the absence of detailed informa- 
tion, one typically sets the value of this common sticking probability to unity, thereby assuming perfect 
electron transmission at the interface. This represents a maximally symmetric situation that is amenable 
to exact mathematical analysis. 

With the simplifications stated above having been made, the ensemble average of a typical squared 
S-matrix clement |5afcp is independent of the channel indices a, b provided a ^ b. Furthermore, in the case 
of the Gaussian Orthogonal Ensemble (GOE), the presence of elastic enhancement (or equivalently — in 



mesoscopic physics — coherent back-scattering) ensures that jS'aal^ = 2|S'abp for all a ^ b. Application 
of the supersymmetry method |^] to problems with GOE symmetry in general yields analytical results 



for \Sab\^ in terms of superintegrals over a 16-dimensional coset-manifold of 8 x 8 supermatrices Q. The 
superintegration ranges over eight commuting (real) degrees of freedom and eight anti-commuting (Grass- 
mann) variables. These superintegrals can normally be reduced further to characteristic triple (ordinary) 
integrals over one compact parameter and two non-compact (semi-infinite) ones. These parameters can 
be identified (modulo coordinate transformations) with the distinct eigenvalues of the supermatrices Q. 
For the problem at hand, one obtains the expression given in Eq. (C.13) of Ref. namely 



Sab\' = {hi , (1.2) 



provided a ^ b, where Im denotes the integral 



(mi - Ai)^(M2 - 
We should note that the change of variables 



2 1 1 



(1.3) 



fi^l- H , ni-l , Ai2^M2-l (1-4) 

has been made relative to Ref. |^. 

It will be advantageous to make the further change of integration variables 

111 . ^ 

M = T, Mi = :r' M2 = — , (1-5) 

so that we can rewrite Im as 

h,-j^ dtj^ dt,J^ "'^[(l-iO(l-t2)]V2-^^-(^-tl)^(^-t2)^-^*"'^^' ^'-'^ 
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where we have made use of the fact that the integrand of Im is symmetric in the variables /ii and ^2- 



By appeaUng to the general principles of (i) unitarity of the S-matrix and of (ii) elastic enhancement 
in the presence of orthogonal symmetry, one can argue that the result 



= 1/(A/+1) (1.7) 

must hold, assuming Tc = 1 for all c = 1,2, . . . , Af , as we have done 1^. Evaluating the integrals in 
Eq. ( |l.3| ) directly is quite another matter, in spite of the simplicity of the final result. Let us also note 
here that the ensuing expression for the mean dimcnsionlcss conductance 'g under the present assumptions 
is 



g{M) = -Y\SabV = ^/m. (1.8) 



The integral in Eq. (1.3) has been evaluated explicitly in Appendix C of Ref. ^ after a lengthy and 



convoluted calculation using a method that will not readily generalize upon the introduction of additional 
complications. One such complication that is of physical interest is the incorporation of a parameter s 
which measures the breaking of orthogonal symmetry due to the presence of a magnetic field and which 
drives the transition towards behaviour described by the Gaussian Unitary Ensemble. If one were to 



imagine expanding the corresponding result for \Sab\'^ in powers or inverse powers of s, in order to study 
the limiting cases of small amd large magnetic fields, respectively, then the coeffients in such an expansion 
would also be triple integrals of the same general nature, which one could expect to evaluate analytically 
as algebraic expressions in the number of channels M . 

Such terms were evaluated using laborious ad hoc procedures in Appendix D of Ref. ^. It would be 
very useful to have a general systematic method for performing these calculations, as well as those for the 
higher dimensional integrals that would arise from considering higher order correlators (like four-point 
functions) in the pure GOE plus external channels problem. 



1.2 Dotsenko-Fateev Integrals 

The integral in Eq. ( |1.6D actually bears a close relationship with a class of integrals introduced first in 
the context of four-point correlation functions in two-dimensional conformal field theories by Dotsenko & 
Fateev . The Selberg integral, which is already known to be intimately connected with the Calogero- 
Sutherland model and with random-matrix theory in the guise of the pure GOE problem without coupling 
to external channels , is in fact a special case of a Dotsenko-Fateev (DF) integral. 

After appropriate variable transformations chosen to cast it into a form suitable for the present 
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application, the general DF integral can be written as 

,{{t,},{rj};a,l3,p)\ , (1.9) 



Jnm{a, 13; p) = — - Y[ / dU Yl / dTj |/„ 



where 



= n^?'(^« - ifU^ri^-^jf - u- r '^2 ■ (1-10) 

t=l f=l n^=lIlJ = lit^ - ^jV 

The constraints on the parameters that ensure solvability of the integral are given by 

a' = -p'a, (3' = -p(3, p = p^^ . (1.11) 

This is the form that appears in Rcf. |l^ where the result of its evaluation can be found. According to 
Eq.(3.8) of Ref. |l^ the value of this integral can be expressed as 

j_(a,/3,p) = p YY^^^YY—^ 

V{\- (3/p + l/p)V {-l + 2m + a/p + /3/p -{n-l + l)lp) 



m — 1 



J=0 



r(a/p - ijp) 

TT r(l - n + g + jp)r(l - n + /3 + jp) 

AJ- r(2-n + a + /3+(m-l+j» ' ^ ' ' 



having implemented the conditions on the parameters a' ,(3' ,p' given in Eq.(l.ll) 



As it turns out, a knowledge of the dependence of the DF integral J7i2(q!,/3; p) on arbitrary values 



of its parameters a,(3,p is required in order to evaluate the particular integral expression (1.6) for Im- 
Specifically, we shall need to consider the more general class of integrals defined by 

J[f] = / dte\t~lf dhf^il-hf dht-^il-hf 

Ji Jo Ja 



\tl~t2\^ 



:f{tMM). (1.13) 



{t-hf{t-h)- 
and we shall require that 

a'^-'i. ..'^-I^ (1,14, 

Then, for /(t, ti, t2) = 1, the integral J[f] belongs to the class of integrals discussed and evaluated analyt- 
ically by Dotsenko & Fateev, i.e., those encompassed by Eqs. (|1.9D-( |lTi] ). We have J[l] = 2j7i2(ck, /3; p), 



and we satisfy the solvability conditions of Eq. (1.11) with the choice 



,'-1- ,-1- (1.15) 

We note however, that the parameter p does not actually enter our expression because we are considering 
the special case where there exists only a single non-compact integration variable (viz. t). 
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What we have encountered in the present appHcation is a moment of the DF integral. We must 
choose 

f{tMM) = t-h , (1.16) 

and 

a = Af/2 , a' = 

/3--1/2, /3' = 1. (1.17) 

These latter relations imply that ultimately A = 1 must be considered. The arbitrariness in A at this 
stage is required to ensure the convergence of all intermediate steps. Thus, in the framework of Dotsenko 
& Fateev, we have p — 1/2, p' — 2, which are values characteristic of problems with GOE symmetry. 

As mentioned above, one is interested in computing moments of DF integrals, such as the one defined 
by Eqs. ( 1.13D and (1.14), in general since they will arise when considering perturbed GOE's and making 



asymptotic expansions of the exact integral result in powers of some small parameter. 



1.3 The Method of Aomoto 



We shall approach the problem at hand by method invented by Aomoto to solve the Selberg integral. 
For this purpose, we introduce the function 

1/3 1^1 ~ ^2!^ 



Dx{t, h, t2) = r{t- ly ihhrni - ii)(i - h)] 



Then 



dt 



a' P' 



(1.18) 



(1.19) 



t t-l t-ti t-t2_ 

Using the symmetry of the function D\{t,ti,t2) in the variables ti and t2, and the fact that it vanishes 
sufficiently rapidly at the boundaries of the t-integration, viz, t = 1 and 00, we see that 



/OO /• 1 /• 1 o 

dt dh dt2j^{tDx) 



dt [ dti [ dt2 
1 Jo ^0 



. 1 



(l + a' + /3')l + /3': — - - 4 



t-l t~t^ 



which implies the identity 

0^ {1 + a' + (3')J[1]+ f3'J 
Likewise, after some algebraic manipulation, we obtain 

/OO /• 1 /• 1 o 

dt dti dt2-^^{t'Dx) 



t - 1 



4 J 



t - ii 



(1.20) 



(1.21) 



00 /• 1 /• 1 

dt I dti / dt2 
1 Jo Jo 



(2 + a' + p')t + P'l + p 



. 1 



t-1 t-ti 



(1.22) 
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which leads to a second identity 

= (2 + a' + (3')J[t] + fi'J[l] + /3'J 



1 



t- 1 



-4J 



t - ti 



Now, subtracting Eq. (1.21) from Eq. (1.22) serves to eliminate the terms involving J[l/(t 
yields the relationship 



= (2 + a' + (3')J[t] - (1 + a')J[l] + 4J 



t - ti 



4J 



t - ti 



Next, we note that 



a P 



ti 1 — ti t — ti t\ ~ t2 



Thus, upon partial integration 

= 



/ dt f dti f dt2 ^itiDx) 
Ji Jo Jo Oti 



dt [ dti [ dt2 
1 Jo Jo 



where we have made use of the fact that (owing to the ti ^ ^2 symmetry) 



OO nl pi 

dt / dti / dt2 

1 Jo Jo - *2 



i / dt I dti I dt2 

" 1 Jo Jo 

OO /•! /•! 

i / dt I dti J dt2Dx 



ti 



tl — t2 t2 — ti 



This yields the identity 



Also, 



0= (l + iA + a + /?)J[l]-/3J 



/.OO /.I /.I Q 

= j dt j dti j dt2 of^itlDx) 



l-tl 



2J 



t - tl 



dt / dti / dt2 
1 Jo Jo 



where we have now used the fact that 



(2 + a)ti - /3-\- + 2-\- + Xti 
i — tl t — tl 



Dx 



dt / dti / dt2 — Dx = \ I dt I dti I dt 



„ 2 / -"2 

1 JO Jo '-1 '-2 Jl JO JO 



''2 



tl — t2 t2 — tl 



Dx 



dt f dti f dt2tiDx 
1 Jo Jo 



After some further algebraic simplification, we arrive at the identity 



0=(2 + A + a + l3)J[ti] + f3J[l] - (iJ 
Subtracting Eq. ( |l.28| ) from Eq. ( |l.3l| ), we find 



1 



1-ti 



2J 



t - tl 



= (2 + A + a + p)J[ti\ - (1 + iA + a) J[l] - 2J 



t - tl 



2J 



i - tl 



which agam has served to ehminate the terms involving J[l/{t— 1)]. To ehminate the terms involving 
J[ti/{t — ti)] and J[t\/{t — ti)], let us now consider the addition of Eq. (1.32) to one half times Eq. (1.24), 
which yields the relation 

= (i(a' + P') - 1) J[t] + (A + a + l3)J[ti] - (i(A - 1) + a + \a') J[l] . (1.33) 

It is more convenient to re-express this relationship as 

= (i(a' + /3') - 1) J[f - 1] + (A + a + (3)J[h - 1] + {\{\ - 1) + /3 + \I3') J[l] . (1.34) 



Next, we must set a' — — 2a/A, j3' = — 2/3/A, which leads to 

Z{a + /5 + A) 



-^[1], (1.35) 



Af - 1 + 2A 

having performed the substitutions a — M/2^ (3 — —1/2 in the last line. Therefore, assuming that the 
integrals implicit in J[t — 1] and J[ti — 1] converge when A ^ 1+, we see that 

Jx=i{t - ti] = ^lim (1 - A)V;,[1] . (1.36) 

We can perform the DF integral J\[l] for A in the neighbourhood of unity by appealing to the results 



of Ref. If we set m = 2, n = 1, p = A/2, then Eq. (3.8) of Ref. |12| (or equivalently Eq. ( |1.12D above) 
implies that 



Accordingly, 



A^i+ (A -1)2 



Mt-h] = (1.38) 



in which case (recalling Eq. (1.2)), 



\Sa,\' = (1-39) 
for a 7^ 6, in agreement with the discussion of Ref. |9[ 

It is important to appreciate that because of the singularity in J\[l] at A = 1, the calculation could 
not have been performed without a knowledge of the DF integral in a general neighbourhood about 
A = 1. Let us close the discussion with a remark on the physical significance of the parameters a, (3, p. 
The value p — ^ is, characteristic of all problems exhibiting an underlying GOE symmetry. In pure GOE 
problems (with no external channels), one typically deals with a — — 1/2,/3 = —1/2. For example, the 
density-density correlator p{Ei)p{E2) is given by Eq. ( |2.10| ) below with x = Ei - E2, t = 0, A = 1/2, 



p = I, q = 2 and l/po equal to the mean level spacing. After some simple changes of variables, the 
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integrand in Eq. (2.10) is seen to be of tlie DF type (1.13) witli a = —1/2, /? = —1/2 and p = |. Thus we 
see that the introduction of M external channels promotes a to a general value. It is an open question 
to determine what sort of additional degrees of freedom should be incorporated in the random-matrix 
model in order to move /? away from its default GOE value. 



2 Calogero-Sutherland Model 



The Calogero-Sutherland model (CSM) describes the quantum mechanics of particles on a line interacting 
through a 1 /r'^ pairwise potential. It has been shown that this model is equivalent to a one-dimensional 
ideal gas of non-interacting anyons, and its excitations are known to exhibit fractional statistics. For 
arbitrary rational values of the coupling X — p/q, p, q £ an exact analytical expression has been given 
by Ha jl^ , in terms of a (p -t- g)-dimensional integral, for the dynamical ground-state density-density 
correlation function (0|/3(a;, t)p(0, 0)|0). The physical significance of the integers p,q is that an excited 
state of the CSM involving p elementary quasi-particle excitations is always accompanied by q quasi- 
hole excitations. The CSM is closely related to the circular ensembles of random matrix theory, with 
the ground-state wavefunctions of the CSM reproducing the eigenvalue distributions of random matrix 
ensembles corresponding to orthogonal, unitary and symplectic symmetry at couplings of A = 1/2, 1,2, 
respectively. Also, after a Wick rotation in the time variable t, the dynamical density-density correlation 
at these couplings coincides with the dynamical density-density correlation for the parameter-dependent 
GOE, CUE and GSE, respectively |l|. It is not too surprising then, to discover that the DF integral 
has a role to play in the calculation of physical quantities in the CSM. 

Starting with the general form of the DF integral, we shall consider a limiting form obtained as the 
result of a certain scaling. This will lead us to an exact identity which allows one to directly compute the 
value of the multi-dimensional integral that represents the overall coefficient in the asymptotic behaviour 
of Ha's density-density correlation function in the limit of large spatial and temporal separations. 



2.1 An Integral Identity 



Suppose that in the DF integral, as given by Eqs. (1^) and ( 1.10 ), we set ti — exTp{si/a} and 



Tj = exp{— ^j/a} after eliminating the parameters a',P',p' with the aid of Eq. (1.11), and then con- 
sider the limit a — > oo. In this limit, we can write 
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i=l j=l 



n:=in;ii«-^""(^.+o)' 

Noting also that 



we see that Eq. (|l.9| ) becomes 

roc ^ PCO 



(2.1) 



1 /"^ /■ 1 

dU ^ - ds, , dTj ^ - I dij , (2.2) 



(2.3) 



where 

CnmiP, p) = m + n + mlS — nP/p + m{m — l)p + n(n — 1)/ p — 2mn . (2-4) 



On the other hand, using the relation 

r(a + x) 



r(a;) 
we see that 

T{l-n + a + jp) r(a) 



X- , (2.5) 



r{2-n + a + (3+{m-l+j)p) a^oo T{a + 1 + (3 + {m - l)p) 

0,^00 Q,l+/3+(m-l)p ' (2-^) 

and 

r(-l + 2m + g/p + /3/p - (n - 1 + £)/p) _ T{a/p - 1 + 2m + - (?i - l)/p) 

T{a/p — i/p) Q^oo r(a/p) 

/ps l-2m-/3/p+(n-l)/p 

^ . (2.7) 



It follows from Eq. (1.12) that 



c„™(/3.p)„(l-/3/p+(n-l)/p)n 7T ^i^/p) TT ^iJP " 

lr(i/p)li r(p) 

n — 1 m — l 

Yl r(l - p/p + II p) [] r(l - n + /? + jp) . (2.8) 



^=0 i=o 



After equating Eq. (2.3) with Eq. ( |2.8| ) in the limit a ^ cxd, we finally obtain the integral identity 



1 " f-OO "'■ f-OO 

m!n! J^JL ± ± j 



n—1 m — l 

' ■ n ^) n n - + ti,) n ni - - + ,3 + 



c„„(/5,p) TT r(^/p) TT r(jp - n) 

Air(i/p)AJ- r(p) 

(2.9) 
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having also introduced the common scahng factor z and used it in making the change of integration 
variables sij p — zr^, = zr\j. 



2.2 Ha's Result 



The exact result of Ha jl^ for the dynamical ground-state density-density correlator reads 



(0|p(a;,t)p(0,0)|0) 



1=1 



J]^ + A)]^ ^ n t-^^-J^^ ^ y^)]^^^ ^cosQa;exp{-i£'i} , 



(2.10) 



i=l 



where 



Q 



1 



c 



(2^Po)' ( 5Z + A) + ^ Ayj(l -yj) \ , 

A2p(9-i)r2(p) r«(A)rp(i/A) 



(2.11) 



2^V<z! nti r2(p - A(* - 1)) n^^i r2(i - (j - i)/A) ' 

and p, q are related by p/ q = A. One should note that, in the formulae given above, the symbol po denotes 
the average density of states for the system po = ^ I L [N being the total number of particles and L the 
one-dimensional volume of the system) , and should not be confused with the parameter of similar name 
p in the DF integral. 

In the region of large x and t, the integral above is saturated by neighbourhoods of all points satisfying 
Xi = Q for all i and ?/j = or 1. However, any point with yj = 1 for some j = 1,2, ... ,p will give rise to 
an oscillating contribution. Thus, in order to calculate the leading-order non-oscillatory term in the limit 
of large x,t, we need to expand all the variables in the integrand about the end-points Xi,yj = 0. This 
yields 

1 1 



(0|p(x,t)p(0,0)|0) 



U27rpo)^CA{q,p,X) 



{27Tipox — i(27rpo)^At)^ {2nipox + i(27rpo)^Ai)^_ 

(2.12) 



where 



A{q,p, A) 



*? poo P pOO 

n/ dx.xf~^e~^^l[ dy,y 



l/\-l^^y, ( _Q_ 

2ttpq 



nLim=i(^.+Ay,)2 



(2.13) 
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In carrying out the final steps in deriving Eq. (2.12), the yj integration intervals have been extended 
to infinity and then changes of integration variables have been made in order to remove the x and t 
dependent factors from the exponentials. Additionally, since the exponentials were initially imaginary, a 
Wick rotation is also assumed to have been performed without obstruction from significant singularities 
in the complex Xi and yj planes. The fact the X = p/q has been used too. 



2.3 Asymptotic Limit 



The task now is to compute the multiple integral A{q,p, A) of Eq. ( 2.15 ). It is clear that differentiating 
the integral identity of Eq. ( |2.9|) twice with respect to z and subsequently setting z = 1 will yield an 
integral of the required type. We must choose f3 = p— 1, m = q, n— p and p = {p + e)/q, with the limit 
e — * to be taken later. Then, 



(2.14) 



■A{q,p,X ^p/q) -^Ipq{z;p-~ l,p) 

having let e ^ on the RHS. Now, to leading order in e, we find 

^-cp,(/3,p) _ ^-e'iq/p) (2.15) 
e->0 

Thus, 

recalling that p = (p + e)/q. After differentiating twice with respect to z, setting z — 1 and taking the 
limit e ^ 0, we obtain the evaluation of the integral which occurs in the asymptotic expansion, i.e.. 

We note that the integral identity (|2.9[) could not be used directly with the final desired choice of 



parameters p = X = p/q because, as seen from Eq. (2T£ ), this case is singular. Only after the twofold 



differentiation can the final parameter values be set. Hence, as in the previous calculation, a knowledge 
of the DF-integral based identity ( ^.9| ) is required in an arbitrary neighbourhood of the relevant point in 
the parameter space. 



Let us now express the normalization constant C in Eq. ( ^.11 ) in a form similar to the product of 
gamma functions appearing above. We have 

1 x-^p{Q-i) 1 

nLiT^b-M^-i)) ^ tHp) ■nL2r^(-A(*-i)) 

1 



ni=2 Hp - - - 1)) Hp - 2)/A - (^ - i)) • ■ • (-(» - i))! ' 
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^_2p(g-l) ^ 

rHp)THqy Ut2^H-H^ i)) }J 



r((j-i)/A-(g-i)) 



r((j ~ i)/A) 



But, 



nr^((j-i)/A-(g-i)) = nr'(i-(j-i)A) 

1^1 ^1 

r^5r2((j-i)/A) = nf^(-7A)' 



and 



9-1 



nr^-A(z-i)) = nr'(*A-p), 

having made extensive use of the relation A = p/ q- Therefore, 
1 A-2p(9-i) 1 



n 



r2(i-(j-i)/A) 



Hence, 



nLir2(p-A(*-i)) r2(p) n?=ir2(*A-p)/i r^o/A) 

r«(A)rp(i/A) 



2^V9! m=l^2(J7A)^?^l^2(^A-p) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



Substitution of this formula for C into the asymptotic expansion above (i.e. Eq. (2.12)), together 



with the evaluation of the multiple integral y^((7,p. A) in Eq. (2.17), gives 

1 1 



i|p(x,Op(0,0)|0) 



-P1 

x,t — >^oo A 



(27rpoa; - (27rpo)2At)2 ' {2TTpax + {2Trpo)'^Xt)'^ 



(2.23) 



which is the required result, i.e. Eq. (7.10) of Ref. |l^. In that paper, the relationship between the 



normalization constant C and the integral A{q,p, A), which are multiplied together in Eq. (2.12) to yield 
the simple factor CA = (27r2A)^^, had to be deduced by indirect means. Here, with the aid of the DP 
integral, we have constructed an explicit derivation of this result. 
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